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INTRODUCTION 

It  is  generally  accepted  that  comparison  of  pre-  versus  post-gadolinium  (gadopentetate 
dimeglumine)  enhanced  gradient  recalled  echo  MRI  offers  exceptionally  high  sensitivity  for 
the  detection  of  breast  cancers  [1-6].  The  specificity  of  gadolinium-enhanced  MR 
mammography  to  distinguish  benign  from  malignant  tissues,  however,  has  been  limited  by 
the  fact  Siat  some  normal  and  benign  masses  also  exhibit  contrast  enhancement  [3,5,6]. 

The  "temporal  rate"  of  gadolinium  enhancement  has  also  generated  considerable  interest  as 
a  potential  complimentary  discriminator  of  benign  and  malignant  breast  neoplasm  [7-10]. 
The  underlying  mechanism  of  contrast  enhancement  rate  as  a  diagnostic  parameter  relates  to 
tissue  vascularity  which  has  been  shown  to  be  a  key  correlate  in  die  growth  and  metastatic 
potential  of  breast  malignancies.  Non-invasive  methods  to  detect  and/or  discriminate 
between  normal,  benign,  pre-malignant  and  malignant  tissues  based  on  vascularity  hold 
inherent  diagnostic  promise  by  virtue  of  this  link  to  tumor  angiogenesis. 

Acquisition  and  processing  methods  to  quantify  enhancement  rate,  or  temporal  properties  in 
general,  have  been  widely  variable  across  studies  [7-10].  In  the  assessment  of  temporal 
enhancement  features,  critical  parameters  include  temporal  resolution,  volume  of  tissue 
coverage,  and  die  method  to  quantify  signal  evolution.  Information  regarding  kinetics  will 
be  obscured  if  the  temporal  resolution  is  inadequate  to  sample  rapid  signal  changes  or 
enhancement  curve  shape  details.  Problems  also  arise  if  the  dynamically  imaged  volume 
only  partially  includes  an  enhancing  lesion  or  does  not  include  all  lesions  in  the  breast. 
These  issues  arise,  in  part,  because  of  the  inevitable  trade-off  in  standard  MRI  acquisition 
techniques  between  volume  coverage,  spatial  resolution,  and  temporal  resolution.  Central 
to  this  project  are  refinement  and  clinical  application  of  new  methods  designed  to  partially 
defeat  this  tradeoff  in  providing  full  volume  coverage  of  both  breasts  at  good  spatial  and 
temporal  resolution  for  detailed  quantification  of  dynamic  contrast  enhancement  [10].  The 
approach,  referred  to  as  "3D  keyhole"  involves  specialized  software  for  data  acquisition 
(i.e.  MRI  pulse  sequence)  and  reconstruction.  This  methodology  allows  dynamic  scans 
through  32  sections  of  both  breasts  to  be  acquired  at  13second  temporal  resolution 
(compared  to  commonly  used  ^minute  resolution  in  volume  scans).  A  major  practical 
advantage  of  this  technique  is  the  ability  to  retrospectively  select  and  analyze  multiple 
arbitrary  tissues/lesions  within  either  breast  that  would  otherwise  not  identified  prior  to 
contrast  administration.  Extensive  data  processing  tools  have  been  developed  to  visualize 
dynamic  contrast  changes  and  reduce  temporal  enhancement  curves  of  lesions  to 
numerically  fitted  time  and  amplitude  constants.  This  methodology  is  being  applied  in  a 
prospective  study  of  100  women  in  whom  a  breast  abnormality,  detected  by  conventional 
means,  will  be  characterized  by  surgical  biopsy. 
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METHODS 

Technical 

Breast  MR  exams  are  performed  on  either  of  our  two  1.5T  General  Electric  Signa  systems 
using  a  commercial  bilateral  phased-array  breast  coil  with  insert  pads  for  mild  A/P 
compression.  After  a  quick  locator  scan,  conventional  axial  T2-wt  fast-spin-echo  is 
performed  to  identify  cysts.  Axial  pre-  and  post-gadolinium  Tl-wt  spin-echo  bracket  the 
"dynamic  scan"  described  below. 

Dynamic  scans  consist  of  three  segments  all  within  one  series  prescription: 

a)  Pre-Gad  Reference:  A  coronal  3D  rf-spoiled  gradient  recalled  echo  (3D  SPGR)  is 
prescribed  to  encompass  both  breasts.  Parameters  are:  TR=13ms;  TE=4.5ms; 
flip=40°;  NEX=4;  typical  FOV=30cm  (R/L),  15cm  (S/I),  3mm  slice  thk;  matrix=256 
(R/L),  128  (S/I),  32  slices  (A/P);  acquisition  time=3'30".  Raw  time-domain  data  is 
automatically  written  to  disk. 

b)  Dynamic-Gad:  A  flag  is  set  for  acquisition  of  central  256  (kR/L)  x  32  (ks/i)  x  32  (kA/p) 
k-space  (3D  keyhole)  data  in  20  serial  passes  at  13seconds  each  pass  (NEX=1).  A  full 
resolution  matrix  "21st  pass"  is  automatically  concatenated  at  the  end  for  an  acquisition 
time  of  519".  A  lOsec  bolus  injection  of  O.lmM/kg  of  gadolinium-DTPA  (Magnevist) 
is  administered  with  4'45"  remaining  in  the  scan.  Rawdata  is  automatically  written  to 
disk. 

c)  Post-Gad  Scan:  A  flag  is  set  to  repeat  scan  with  parameters  in  listen  in  a). 

The  total  MR  exam  time  is  typically  40-45  minutes. 

Raw  data  is  networked  to  a  SUN  SPARC10  for  offline  3D  keyhole  processing.  This 
involves  some  corrections  and  insertion  of  each  keyhole  dataset  into  the  full  matrix 
reference  dataset,  then  3D  Fourier  transform  to  yield  32  slices  at  each  of  20  timepoints. 

The  21st  pass  dataset  is  used  to  correct  an  enhancement  amplitude  error  that  results  from 
partial  k-space  acquisition  [10].  Keyhole  reconstructed  images  and  "reference-subtracted" 
images  are  viewed  cinegraphically.  Enhancing  lesions  are  best  perceived  in  this  format 
which  guides  ROI  definition.  Lesion  enhancement  time-course  curves  are  then 
automatically  fit  to  mono-  and  double-exponential  saturation  models  such  that  temporal 
features  are  reduced  to  enhancement  time-constants  [10].  Additional  software  has  been 
written  to  cinegraphically  display  maximum-intensity-projection  (MIP)  of  subtracted 
images  through  the  full  3D  volume  and  display  heterogeneity  of  fitted  constants  within  the 
defined  ROI. 

Clinical 

Recruited  patients  are  referred  from  the  Breast  Imaging  Section  of  the  Breast  Care  Center  at 
the  University  of  Michigan  Hospital.  Entrance  criteria  includes  any  patient  with  a  breast 
abnormality  detected  by  mammography,  physical  examination,  or  ultrasound  which  will  be 
biopsied.  Prior  to  MRI  study,  all  patients  have  a  routine  clinical  evaluation  which  includes, 
as  a  minimum,  mammography  and  breast  physical  examination.  Twenty-five  patients  were 
to  be  studied  in  each  of  4years.  Twenty-seven  patients  have  been  studied  to  date  in  Yearl. 

Three  radiologists  independently  review  the  data  from  the  following  "modalities"  in  the 
order  of:  (a)  conventional  mammography  images  (Xray,  ultrasound)  and  reports;  (b) 
conventional  MRI  including  pre-  and  post-gadolinium;  and  (c)  subtracted  MR  images  (pre¬ 
gad  subtracted  from  post-gad).  Readers  grade  lesions  on  five  point  scales  of  suspicion 
(benign  to  malignant)  and  visibility  (obvious  to  subtle). 

The  temporal  enhancement  properties  of  tissues/masses  are  quantified  by  numerical  fit  to 
the  21  points  of  the  enhancement  curve  to  yield  mono-  and  double-exponential  saturation 
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time  constants  and  enhancement  amplitudes.  Time-course  data,  fitted  curves  and  related  fit 
parameters  are  stored  on  disk  as  a  permanent  ROI  record. 

RESULTS 

Determination  of  the  diagnostic  power  of  temporal  contrast  enhancement  properties  as 
quantified  by  time-constants  is  central  to  the  goal  of  this  research.  We  have  combined  data 
from  Year  1  (27  patients)  to  that  of  our  pilot  study  using  similar  methodology  to  increase 
the  study  population  to  41  patients.  A  scatterplot  of  fitted  mono-exponential  time  constant 
versus  tissue  type  is  illustrated  in  the  figure. _ 
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Figure:  Enhancement  time 

, _ constant  vs  tissue  type _ 


Tissue  Type  Key: 


1  =  Inv  Ductal  Ca 

2  =  Neuroendocrine  Ca 

3  =  Inv  Lobular  Ca 

4  =  IntraDuctal  Ca 

5  =  Adj.  to  InvDuctal  Ca 

6  =  Lobular  Ca  InSitu 

7  =  Fibroadenoma 

8  =  Fibrocystic  change 

9  =  Benign  mass 
10=  Nipple  /  Areola 
11=  Normal  parenchyma 


CONCLUSIONS 

Study  hypotheses  relate  to  the  diagnostic  value  of  contrast  enhancement  temporal  features 
in  distinguishing  benign  from  malignant  and  between  tissue/tumor  types.  The  hypotheses 
for  testing  state  that  a  high  rate  of  contrast  change  (i.e.  short  enhancement  time  constant) 
raises  suspicion  of  malignancy.  While  the  data  generated  thus  far  are  limited  by  patient 
numbers,  the  trend  for  malignancies  to  exhibit  rapid  enhancement  relative  to  non-malignant 
tissues  is  apparent.  Diagnosis  based  solely  on  enhancement  temporal  features,  however, 
does  not  guarantee  against  false  positives  or,  more  seriously,  false  negatives  as  indicated  in 
the  figure.  The  prevalence  of  these  misclassifications  require  greater  patient  volume  as 
planned  for  the  remainder  of  the  study. 
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Dynamic  Three-dimensional  Imaging 
with  Partial  K-Space  Sampling:  Initial 
Application  for  Gadolinium-enhanced  Rate 
Characterization  of  Breast  Lesions1 


PURPOSE:  To  evaluate  a  method  to 
monitor  gadolinium  enhancement  pat¬ 
terns  at  magnetic  resonance  (MR)  imag¬ 
ing  with  high  temporal  resolution  and 
full  coverage  through  both  breasts. 

MATERIALS  AND  METHODS:  In  12 
patients  with  13  masses,  including 
nine  carcinoma,  nonenhanced  three- 
dimensional  MR  imaging  was  per¬ 
formed  with  full-matrix  resolution. 

At  dynamic  imaging,  32  serial  passes 
were  made  during  bolus  administra¬ 
tion  of  contrast  material,  and  tempo¬ 
ral  resolution  was  reduced  to  12  sec¬ 
onds  by  collecting  the  central  (low 
spatial  frequency)  32  x  16  or  16  x  16 
phase-encode  views.  Full-matrix  dy¬ 
namic  images  were  reconstructed  by 
complementing  central  phase-encode 
data  with  precontrast  data  from  pe¬ 
ripheral  high-spatial-frequency  views. 

RESULTS:  Results  at  time-course 
analysis  with  a  mono-exponential 
saturation  model  indicated  malignant 
lesions  tend  to  show  rapid  (<60  sec¬ 
onds)  contrast  change  relative  to  be¬ 
nign  masses  and  normal  tissues.  One 
cancer  displayed  an  exceptionally 
slow  contrast  change  (260  seconds). 

CONCLUSION:  The  technical  objec¬ 
tives  of  full  tissue  coverage,  rapid 
temporal  sampling,  and  quantifica¬ 
tion  of  enhancement  curves  are  met 
with  this  method  for  certain  lesions 
(>5  mm  in  largest  diameter). 

Index  terms:  Breast,  MR,  00.12143  •  Breast 
neoplasms,  MR,  00.12143  •  Magnetic  resonance 
(MR),  k-space 
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Early  magnetic  resonance  (MR) 
imaging  demonstrated  only  lim¬ 
ited  success  in  assisting  diagnosis  of 
breast  disease  (1-4).  Initially,  only  tis¬ 
sue-inherent  differences  in  nuclear 
MR  relaxation  times,  T1  and  T2,  were 
the  primary  sources  of  contrast.  Al¬ 
though  these  contrasts  are  relatively 
high,  it  was  observed  that  they  have 
little  value  in  helping  detection  and 
diagnosis  of  breast  cancer  (1,3,5,6). 
Variable  admixtures  of  intense  fat  sig¬ 
nal  and  inadequate  fat  suppression 
have  also  seriously  hindered  image 
interpretation  and  detection  of  breast 
carcinoma. 

More  recently,  gadolinium-en¬ 
hanced  MR  mammography  has  gen¬ 
erated  considerable  interest  (7-10).  It 
is  generally  accepted  that  gadolinium 
enhanced  MR  imaging  offers  high 
sensitivity  for  the  detection  of  breast 
cancers.  The  specificity  of  gadolinium- 
enhanced  MR  mammography  to  dis¬ 
tinguish  benign  from  malignant  tis¬ 
sues,  however,  has  been  limited  by 
the  fact  that  benign  tissues  can  also 
exhibit  contrast  enhancement  (9,11,12). 

The  rate  of  gadolinium  enhance¬ 
ment  has  also  been  studied  as  a  po¬ 
tential  discriminator  of  benign  and 
malignant  breast  neoplasms  (13-15). 

In  their  initial  summary  of  findings  in 
25  patients  who  underwent  dynamic 
gadolinium-enhanced  MR  imaging, 
Kaiser  and  Zeitler  reported  that  all 
malignancies  displayed  rapid  en¬ 
hancement  (ie,  ~  100%  increase  in 
signal  intensity  within  the  first  2  min¬ 
utes  after  injection  of  0.1  mmol/kg 
gadopentetate  dimeglumine),  whereas 
benign  lesions  demonstrated  less  and 
slower  signal  intensity  change  (13). 


These  authors  have  since  updated 
their  dynamic  gadolinium-enhanced 
MR  imaging  study,  in  which  sensitiv¬ 
ity  and  specificity  for  cancer  detection 
was  greater  than  95%  (16).  In  other 
dynamic  gadolinium-enhanced  MR 
imaging  studies,  higher  false-positive 
rates  were  observed;  thus,  there  is  not 
a  consensus  among  investigators  on 
methods  or  diagnostic  value  of  this 
approach  (11,12,17).  Findings  in  these 
and  other  studies  illustrate  that  fibro¬ 
adenomas  and  proliferative  fibrocystic 
change,  in  particular,  are  sources  of 
false-positive  findings.  Breast  malig¬ 
nancies  typically  exhibit  rapid  en¬ 
hancement  while  benign  lesions  usu¬ 
ally,  although  not  necessarily,  exhibit 
slower  contrast  change.  Malignant 
lesions  may  enhance  slowly,  thus  re¬ 
ducing  sensitivity  (12,18). 

Use  of  various  acquisition  and  en¬ 
hancement  rate  quantification  schemes 
may  contribute  to  discrepancies  in 
findings  (13-15).  In  assessment  of 
enhancement  rate,  important  param¬ 
eters  include  temporal  resolution, 
volume  of  tissue  covered,  and  the 
method  used  to  quantify  signal  inten¬ 
sity  change.  In  particular,  important 
information  on  the  rate  of  enhance¬ 
ment  may  be  obscured  if  the  temporal 
resolution  is  inadequate  to  sample 
rapid  contrast  changes.  In  addition, 
problems  will  obviously  arise  if  the 
dynamic  imaging  volume  includes 
an  enhancing  lesion  only  partially 
or  does  not  include  all  lesions  in  the 
breast.  These  issues  arise  in  part  be¬ 
cause  there  is  an  inevitable  trade-off 
in  standard  MR  imaging  between  im¬ 
aged  volume  and  temporal  resolution. 

MR  imaging  with  methods  such  as 
high  spatial  resolution,  three-dimen¬ 
sional  (3D)  volumes,  suppression  of 
fat  signal,  and  magnetization  transfer 
contrast  have  been  studied  in  an  at- 


Abbreviations:  ROI  =  region  of  interest,  3D  = 
three-dimensional. 
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tempt  to  produce  high-quality  ana¬ 
tomic  images  (9,19).  These  3D  acquisi¬ 
tions  typically  require  several  minutes; 
thus,  dynamic  properties  are  not  readily 
observed.  Conversely,  to  improve 
temporal  sampling  of  enhancement 
curves,  others  have  performed  two- 
dimensional  multisection  MR  imaging 
with  1-minute  temporal  resolution 
(13)  or  of  one  judiciously  selected  sec¬ 
tion  with  finer  temporal  (15-second) 
resolution  (14).  Unfortunately,  selection 
of  prospective  section(s)  for  dynamic 
MR  imaging  by  means  of  correlation 
with  findings  at  x-ray  mammography, 
ultrasonography  (US),  or  physical  ex¬ 
amination  can  be  extremely  difficult 
in  the  amorphous  anatomy  of  a  vari¬ 
ably  compressed  breast,  except  in 
cases  of  large  tumors.  Moreover,  the 
nonspecific  features  of  some  lesions 
on  precontrast  MR  images  may  fur¬ 
ther  hinder  prospective  selection  of 
appropriate  section(s)  for  dynamic 
imaging. 

The  central  objective  of  this  study 
was  to  obviate  the  trade-off  between 
tissue  coverage  and  temporal  resolu¬ 
tion  seen  in  other  studies  of  dynamic 
gadolinium-enhanced  MR  imaging  of 
the  breast.  The  specific  technical  ob¬ 
jectives  were  development  and  vali¬ 
dation  of  a  pulse  sequence  that  would 

(a)  monitor  gadolinium  enhancement 
on  a  short  time  scale  ( ~  12  seconds), 

(b)  image  the  full  volume  of  both 
breasts,  (c)  produce  anatomic  images 
with  sufficient  spatial  resolution  and 
contrast  to  allow  correlation  of  find¬ 
ings  with  those  on  images  obtained 
with  other  modalities,  and  (d)  quan¬ 
tify  parameters  related  to  enhance¬ 
ment  rate. 

MATERIALS  AND  METHODS 
Patients 

Twelve  female  patients,  aged  30-62 
years  (mean,  48  years),  were  recruited  for 
this  study.  Eleven  patients  had  one  or 
multiple  masses  identified  at  mammogra¬ 
phy  and/or  physical  examination.  Patho¬ 
logic  confirmation  was  obtained  in  11  of  13 
masses,  including  nine  invasive  carcinoma 
tumors  (largest  diameter  mean,  22  mm; 
range,  12-50  mm)  in  seven  patients  and  in 
the  one  benign  mass  in  each  of  two  pa¬ 
tients  (fibrosis  with  largest  diameter  of  30 
mm  and  a  palpable  benign  mass  with  larg¬ 
est  diameter  of  10  mm).  Of  the  remaining 
two  of  13  masses,  one  was  a  palpable  mass 
with  a  largest  diameter  of  20  mm  that  re¬ 
solved  with  time  in  one  patient  and  was 
therefore  presumed  to  be  benign  and  the 
other  was  a  mass  with  a  largest  diameter 
of  17  mm  that  was  radiographically  consis¬ 
tent  with  a  fibroadenoma  and  was  stable 
for  7  years  in  another  patient.  Cytologic 
findings  were  also  consistent  with  a  fibro¬ 


adenoma.  One  patient  had  no  detectable 
mass  and  was  presumed  to  have  no  breast 
disease.  The  MR  imaging  methods  were 
approved  by  an  institutional  review 
board,  and  informed  written  consent  was 
obtained  from  all  participants. 


Data  Acquisition  for  3D  Keyhole 
MR  Imaging 

Patients  underwent  imaging  on  a  Signa 
imager  (GE  Medical  Systems,  Milwaukee, 
Wis).  Plane  orientation  and  field  of  view 
for  the  dynamic  series  were  chosen  to  al¬ 
low  efficient  imaging  of  the  breasts  with 
high  spatial  resolution  and  minimal  image 
wraparound  and  involvement  of  the 
heart.  Coronal  imaging  was  performed 
through  the  suspended  breasts  with  use 
of  a  rectangular  field  of  view  with  matrix 
dimensions  of  256  (right  to  left)  x  128  (su¬ 
perior  to  inferior)  pixels  as  the  preferred 
geometry.  In  eight  patients  imaging  was 
performed  with  a  dedicated  breast  coil. 

In  four  other  patients  imaging  was  per¬ 
formed  with  the  patient  lying  prone  on 
a  surface  coil,  as  a  result  of  body  habitus, 
and  axial  dynamic  imaging  was  necessary 
in  two  of  these  patients.  The  remainder  of 
the  dynamic  studies  were  performed  with 
coronal  views. 

The  patients  initially  underwent  local¬ 
ization  with  conventional  nonenhanced 
Tl-  and  T2-weighted,  spin-echo  se¬ 
quences.  Conventional  enhanced  spin- 
echo,  Tl-weighted  imaging  was  also  per¬ 
formed.  A  reduced  k  space  or  "keyhole" 
technique  (20-22),  described  in  detail  be¬ 
low,  was  used  for  dynamic  imaging.  Pre¬ 
contrast  acquisitions,  referred  to  as  refer¬ 
ence  data  sets,  were  performed  with  a  3D 
matrix  of  256  x  128  x  32  sections  and  one 
signal  acquired. 

Reference  imaging  was  followed  imme¬ 
diately  with  dynamic  3D  keyhole  imaging, 
without  changing  the  hardware  settings, 
and  only  the  central  phase-encode  views 
(in  both  phase-encode  directions)  were 
serially  imaged  in  32  passes  during  admin¬ 
istration  of  contrast  material.  In  nine  of 
the  12  patients  studied,  3D  keyhole  dimen¬ 
sions  were  256  x  32  x  16  (ie,  one-eighth 
of  k  space  acquired),  and  in  the  remain¬ 
ing  three  patients  3D  keyhole  dimensions 
were  256  x  16  x  16  (ie,  one-sixteenth  of 
k  space  acquired).  Timing  parameters  in¬ 
cluded  repetition  time  of  24  msec  and 
echo  time  of  4.5  msec  (24/4.5)  for  the 
256  x  32  x  16  keyhole  acquisitions  and 
45/7  for  the  256  x  16  x  16  keyhole  acqui¬ 
sitions.  Temporal  resolution,  defined  as 
the  time  to  complete  each  keyhole  pass, 
was  comparable  in  all  acquisitions  (ie,  12.3 
and  11.5  seconds,  respectively).  Reference 
and  dynamic  keyhole  imaging  were  both 
performed  with  use  of  the  same  radio-fre¬ 
quency-spoiled  sequence  of  gradient-re¬ 
called  acquisition  in  the  steady  state, 
which  included  a  frequency-selective,  fat- 
saturation  pulse.  Other  imaging  param¬ 
eters  included  a  40°  flip  angle,  one  signal 
acquired,  32-bit  digitization,  and  ±32-kHz 
bandwidth.  The  nominal  right-to-left  field 
of  view  of  300  mm  and  section  thickness 


of  3-6  mm  (mean,  4.5  mm)  were  adjusted 
for  patient  size.  The  position  and  thick¬ 
ness  of  3D  slabs  (equal  to  28  sections)  were 
graphically  prescribed  on  orthogonal  im¬ 
ages  to  encompass  the  volume  of  the 
breasts.  Imaging  times  for  reference  and 
serial  3D  keyhole  acquisitions,  respec¬ 
tively,  were  1  minute  39  seconds  and  6 
minutes  33  seconds  with  repetition  time 
of  24  msec,  and  3  minutes  5  seconds  and  6 
minutes  9  seconds  with  repetition  time  of 
4  msec. 

In  all  cases,  0.1  mmol/kg  of  gadopen- 
tetate  dimeglumine  (Magnevist;  Berlex 
Laboratories,  Wayne,  NJ)  was  manually 
injected  as  a  bolus  over  10-15  seconds  in 
an  antecubital  vein.  The  injection  was 
started  5  minutes  30  seconds  before  the 
end  of  dynamic  imaging.  After  dynamic 
imaging,  the  raw  data  were  permanently 
stored  (requiring  1-2  minutes).  Then  a 
full-matrix  series  was  obtained  with  the 
same  imaging  parameters  used  for  the  ref¬ 
erence  data  set  to  enable  comparison  of 
pre-  and  postcontrast  images. 


Data  Processing 

Time-domain  data  from  precontrast  and 
dynamic  acquisitions  were  transferred  to  a 
computer  (SparclO;  Sun  Microsystems, 
Mountain  View,  Calif)  for  off-line  process¬ 
ing.  Custom  software  tools  for  reconstruc¬ 
tion,  display,  and  region-of-interest  (ROI) 
analysis  were  written  for  processing  (Ad¬ 
vanced  Visual  Systems,  Waltham,  Mass). 
Modules  were  also  written  to  determine 
phase  differences  between  the  early  dy¬ 
namic  passes  (ie,  prior  to  changes  in  con¬ 
trast)  and  the  corresponding  central  time- 
domain  data  of  the  full-matrix  reference 
data  set.  Ideally,  these  should  have  identi¬ 
cal  phase.  Slight  patient  repositioning, 
however,  could  result  in  phase  error  and 
artifact  on  reconstructed  keyhole  images. 
Therefore,  phase  differences  were  deter¬ 
mined  and,  when  necessary,  a  single 
phase-correction  setting  was  applied  to 
all  32  temporal  passes  of  each  study. 

After  phase  correction,  each  pass  was 
substituted  for  the  central-time-domain 
core  of  the  full-matrix  reference  data  set 
and  was  reconstructed  into  a  32-section 
image  set  by  means  of  conventional  3D 
Fourier  transform.  This  process  was  re¬ 
peated  for  all  passes  to  produce  1,024 
anatomic  images  (32  sections  x  32  time 
points).  Subtraction  images  were  also  pro¬ 
duced  by  using  an  early  time  point  as  a 
mask.  Anatomic  and  subtraction  images 
versus  time  were  viewed  side  by  side  cine- 
graphically.  Other  display  formats  de¬ 
signed  for  quick  survey  of  data  sets  in¬ 
cluded  an  image  matrix  of  all  sections 
(anatomic  and  subtracted  images)  at  the 
last  time  point  and  cinegraphic  display 
of  maximum  intensity  projections  of  sub¬ 
tracted  images.  The  maximum-intensity- 
projection  images  were  obtained  through 
the  full  tissue  volume  along  orthogonal 
views  and  were  particularly  beneficial  for 
viewing  changes  in  spatial  and  temporal 
contrast  of  the  lesion(s)  relative  to  the 
heart  and  large  vessels. 
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The  time  course  of  signal  change  in  an 
ROI  of  enhancing  tissues  was  fit  to  a 
mono-exponential  saturation  model  with 
use  of  a  Levenberg-Marquardt  least- 
squares  optimization  algorithm  (23).  Fitted 
parameters  include  contrast  onset,  tQ ;  am¬ 
plitude  of  signal  change,  A ;  and  enhance¬ 
ment  time  constant,  t,  used  within  the  fol¬ 
lowing  model: 


AS(0  = 


fA[  1  -  e~{/-'0)/7],  for  t  >  t0 
0,  for  t  <  to. 


(1) 


RESULTS 

Regarding  mass  conspicuity,  con¬ 
ventional  precontrast  spin-echo  Tl- 
and  T2-weighted  images  frequently 
did  not  delineate  location  or  spatial 
extent  of  the  masses.  Visual  compari¬ 
son  of  pre-  versus  postcontrast  images 
(Tl-weighted  spin-echo  and  gradient- 
recalled-echo  images)  was  adequate 
to  identify  nearly  all  (12  of  13)  pal¬ 
pable  and/or  mammographically  de¬ 
tected  masses.  Mass  detection  and 
localization  were  facilitated  by  inspec¬ 
tion  of  anatomic  and  subtracted  cine- 
graphic  images  generated  by  means 
of  the  keyhole  scheme.  Subtraction 
removed  many  image  features  that 
partially  obscured  lesions,  and  cine- 
graphic  display  further  delineated 
lesions  in  terms  of  time  and  intensity 
of  contrast  change.  As  the  cinegraphic 
format  provided  the  greatest  conspi¬ 
cuity  of  enhancing  tissues,  it  was  used 
to  interactively  guide  ROI  definition 
on  the  time-series  images.  One  super¬ 
ficial  palpable  mass  ( « 1  cm  in  largest 
diameter)  displayed  signal  intensity 
enhancement,  although  the  enhance¬ 
ment  was  not  distinct  relative  to  sur¬ 
rounding  tissues.  The  mass  region 
was  identified  as  a  capsule  on  the  skin 
that  was  visible  at  MR  imaging.  This 
mass  was  occult  at  conventional  x-ray 
mammography  and  US  and  was 
shown  to  be  benign  at  pathologic 
analysis. 

Figure  1  illustrates  data  obtained  in 
a  56-year-old  patient  with  pathologi¬ 
cally  proved  invasive  lobular  carci¬ 
noma  in  two  adjacent  sites  separated 
by  approximately  1  cm.  Mammo¬ 
grams  depicted  only  the  larger  mass. 
An  early-time-point  (ie,  before  con¬ 
trast  intensity  change)  keyhole  image 
(Fig  la)  does  not  depict  the  two  ma¬ 
lignant  foci,  which  were  most  appar¬ 
ent  on  the  late-time-point  subtraction 
keyhole  image  (Fig  lb).  Keyhole  ma¬ 
trix  size  applied  in  this  case  was  256  x 
32  x  16.  Several  keyhole  subtraction 
images  from  the  32-image  time  series 
for  this  section  are  shown  in  Figure 
lc.  While  the  relatively  rapid  en¬ 
hancement  period  is  visually  appar¬ 
ent  in  the  time  series,  temporal  con¬ 


trast  patterns  of  various  tissues  are 
best  viewed  graphically,  as  shown  in 
Figure  Id.  Mono-exponential  param¬ 
eters  and  fitted  curves  are  also  shown 
in  this  plot.  The  time  constants  of  the 
lesions  (t  =  52  or  55  seconds)  indicate 
that  the  majority  of  signal  intensity 
change  occurs  in  less  than  1  minute. 

The  number  or  location  of  these 
masses  was  not  readily  apparent  on 
precontrast  images  (Fig  la),  although 
their  region  roughly  correlated  with  a 
suspect  mass  seen  on  the  mammo¬ 
gram. 

An  example  of  a  slowly  enhancing 
lesion  is  Figure  2,  which  was  obtained 
in  a  56-year-old  patient  with  a  prob¬ 
able  fibroadenoma.  Findings  at  MR 
imaging  and  mammography  were 
easily  correlated  for  this  particular 
mass,  which  had  been  sampled  cyto- 
logically  7  years  earlier.  No  change  in 
size  was  seen,  as  noted  in  Materials 
and  Methods.  A  subset  of  the  subtrac¬ 
tion-image  time  series  illustrates 
gradual  enhancement  of  the  lesion 
(Fig  2c).  The  keyhole  matrix  size  ap¬ 
plied  in  this  case  was  also  256  x  32  x 
16.  As  with  the  data  in  Figure  Id,  the 
time-course  curve  is  well  fit  by  an  ex¬ 
ponential  (Fig  2d),  although  for  this 
lesion  the  time  constant  is  substan¬ 
tially  longer  (t  =  119  seconds). 

A  summary  of  fitted  mono-expo¬ 
nential  time  constants  for  a  variety  of 
tissue  types  and  all  patients  is  shown 
in  the  scatter  plot  of  Figure  3.  ROIs 
from  tissues  labeled  "normal"  were 
derived  from  areas  remote  to  the  mass 
that  displayed  sufficient  enhancement 
for  curve  fitting  but  were  otherwise 
normal  in  appearance.  These  were 
typically  regions  of  diffuse  enhance¬ 
ment  from  either  the  ipsilateral  or 
contralateral  breast.  No  more  than 
one  normal-tissue  ROI  was  derived  in 
each  patient.  Two  of  the  five  invasive 
ductal  carcinoma  tumors  were  found 
in  one  patient  (within  one  breast  and 
separated  by  approximately  5  cm), 
and  two  of  the  three  invasive  lobular 
carcinoma  tumors  were  found  in  one 
patient  (Fig  1).  ROIs  of  tissues  adja¬ 
cent  to  invasive  ductal  carcinoma  tu¬ 
mors  were  defined  to  avoid  signal 
blur  from  the  lesion  focus.  As  men¬ 
tioned  in  Materials  and  Methods,  the 
keyhole  matrix  size  was  256  X  16  x  16 
in  three  of  12  patients.  Masses  charac¬ 
terized  in  these  patients  were  three 
invasive  ductal  carcinoma  tumors 
(with  largest  diameters  of  18, 13,  and 
22  mm  and  t  of  14,  20,  and  18  seconds, 
respectively)  and  one  neuroendocrine 
carcinoma  (largest  diameter,  14  mm; 
t,  37  seconds). 

There  was  a  clear  tendency  for  the 
malignant  lesions  to  display  rapid 


contrast  changes  (t  <  60  seconds) 
relative  to  benign  or  normal  tissues. 
However,  one  invasive  lobular  cancer 
displayed  an  exceptionally  long  time 
constant  with  t  of  260  seconds.  Patho¬ 
logic  examination  revealed  the  tumor 
had  a  largest  diameter  of  at  least  1  cm 
(margin  was  positive  at  initial  biopsy), 
gross  size  of  1.2  x  1.0  x  0.8  cm,  no 
angiolymphatic  invasion,  and  no  ex¬ 
tensive  ductal  carcinoma  in  situ;  in 
addition,  17  of  17  lymph  nodes  were 
negative  for  malignancy. 

DISCUSSION 

We  tested  a  technique  with  which 
to  perform  dynamic,  contrast  materi¬ 
al-enhanced  MR  imaging  of  the  en¬ 
tire  volume  of  both  breasts  with  a 
temporal  resolution  of  approximately 
12  seconds.  The  usefulness  of  the 
method  was  explored  in  our  study  of 
12  patients  with  a  total  of  13  breast 
masses.  As  found  in  other  studies,  we 
noted  a  trend  for  malignant  cancers  to 
display  rapid  enhancement  relative  to 
benign  and  normal  tissues,  although 
the  presence  of  one  slowly  enhancing 
cancer  was  observed.  This  would  be 
considered  a  false-negative  case  if  the 
enhancement  time  constant  were 
used  as  the  sole  diagnostic  indicator 
of  malignancy.  The  small  number  of 
patients  in  this  study,  however,  pre¬ 
vents  us  from  drawing  statistically 
meaningful  inferences  regarding  the 
diagnostic  value  of  this  technique. 
Technical  issues  aside,  other  investi¬ 
gators  have  noted  both  false-negative 
and  false-positive  findings  relative  to 
enhancement  rate  (11,12,17,18). 

The  viability  of  this  approach  for 
management  of  breast  disease  neces¬ 
sitates  determination  of  sensitivity 
and  specificity  rates  for  each  applica¬ 
tion.  The  role  that  we  are  currently 
evaluating  is  the  contribution  of  in¬ 
cremental  information  beyond  that 
provided  with  conventional  imaging 
modalities  when  a  suspect  mass  is 
detected.  In  this  capacity,  emphasis  is 
placed  on  characterization  of  previ¬ 
ously  detected  masses  that  are  usually 
more  than  5  mm  in  largest  diameter. 
Lesion  localization  for  eventual  dy¬ 
namic  gadolinium-enhanced  MR  im¬ 
aging  is  hindered  by  nonspecific  le¬ 
sion  properties  in  precontrast  MR 
images  and/or  by  poor  correlation  of 
findings  with  those  of  other  imaging 
modalities,  owing  to  variable  breast 
geometries  and  views.  This  problem  is 
essentially  solved  with  the  3D  key¬ 
hole  approach,  which  has  the  inher¬ 
ent  features  of  complete  coverage 
through  both  breasts  and  high  tempo¬ 
ral  sampling  of  dynamic  contrast 
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Figure  1.  Data  from  a  coronal  image  obtained  through  the  breast  of  a 
56-year-old  patient  with  two  foci  of  invasive  lobular  carcinoma,  (a)  Early- 
time-point  (ie,  precontrast)  keyhole  image  does  not  delineate  the  cancer 
tumors,  (b)  Subtraction  image  from  a  late  time  point  (t  =  355  seconds)  in 
dynamic  imaging  shows  the  two  high-contrast  cancer  tumors  (ROI#l, 
R01#2).  (c)  Subset  of  time  series  for  this  section  illustrates  the  relatively 
rapid  contrast  change  in  the  cancer  lesions.  "  =  seconds,  (d)  Time-course 
analysis  of  regions  RO/#l  and  ROI# 2  with  use  of  mono-exponential 
saturation  model  whereby  both  cancers  are  characterized  by  a  relatively 
short  time  constant  (t  =  52, 55  seconds)  (ie,  a  rapid  contrast  change). 
Images  were  electronically  cropped  to  show  only  the  area  of  interest. 


changes.  These  advantages  would 
also  be  a  benefit  in  other  potential 
applications,  including  use  as  an  alter¬ 
native  breast  screening  modality  in 
problematic  cases  (eg,  dense  breasts), 
as  a  survey  modality  to  search  for  ad¬ 
ditional  malignancies  prior  to  lumpec¬ 
tomy,  as  an  anatomic  guide  to  direct 
biopsies,  and  potentially  as  an  imag¬ 
ing  modality  in  patients  with  sero¬ 
logic  or  genetic  positive  markers  and 
negative  findings  at  mammography. 
Our  encouraging  preliminary  results 
suggest  that  further  investigation 
with  this  technique  is  warranted. 

The  raw-data  domain  of  MR  imag¬ 
ing,  referred  to  as  k  space,  defines  the 
spatial  frequency  distribution  of  an 
object  and  is  a  convenient  format 
with  which  to  describe  various  data 
acquisition  schemes.  In  a  simplistic 
sense,  gross  contrast  features  are  pri¬ 
marily  determined  by  data  near  the 
so-called  center  of  k  space  (ie,  low 
spatial  frequencies),  whereas  features 
with  high  spatial  detail  are  defined  by 
data  from  more  peripheral  k  space. 
These  concepts  are  used  with  great 
success  in  techniques  that  generate 
desirable  contrast  properties  more 
efficiently  than  routine  MR  imaging 
does  (24,25). 

Since  most  changes  in  object  and 
temporal  contrast  are  centralized  in 
k  space,  other  investigators  have  in¬ 
creased  temporal  resolution  in  dy¬ 
namic  gadolinium-enhanced  two- 
dimensional  keyhole  MR  imaging  by 
repeatedly  sampling  only  the  central 


phase-encoded,  ky,  lines  (20).  Images 
with  high  anatomic  quality  and  tem¬ 
poral  resolution  are  produced  by 
combining  dynamic  data  with  previ¬ 
ously  acquired  peripheral  phase-en¬ 
coded  data.  In  application  to  the 
breast,  a  3D  acquisition  is  highly  de¬ 
sirable  for  its  advantages  in  tissue 
coverage,  spatial  resolution,  and  sig- 
nal-to-noise  ratio.  The  second  phase- 
encode  dimension,  kz,  has  an  image¬ 
time  penalty  proportional  to  the 
number  of  sections.  Fortunately,  the 
keyhole  concept  may  also  be  ex¬ 
tended  to  this  dimension  so  rapid  dy¬ 
namic  volume  imaging  is  possible  (21, 
22).  One  version  of  this  sequence  dy¬ 
namically  acquired  only  one-sixteenth 
of  k  space  and  was  applied  in  three  of 
12  patients  in  this  series;  for  the  re¬ 
maining  studies,  one-fourth  of  ky  and 
one-half  of  kz  (ie,  one-eighth  of  k 
space)  were  dynamically  acquired  in 
one-eighth  the  image  time. 


The  main  drawback  to  sampling 
only  a  portion  of  k  space  is  reduced 
spatial  resolution.  This  is  mitigated 
because  dynamically  acquired  key¬ 
hole  data  can  be  complemented  with 
previously  acquired  full-matrix  data 
to  generate  high-resolution  anatomic 
images.  We  emphasize,  however,  that 
temporally  changing  features  remain 
blurred  by  an  amount  dependent  on 
keyhole  matrix  size  and  lesion  spatial 
dimensions.  That  is,  spatial  resolution 
in  contrast-enhancing  structures  will 
be  limited  by  the  keyhole  resolution 
and  not  by  the  full-matrix  resolution. 
The  rationale  for  splicing  keyhole 
data  into  full-matrix  reference  data  is 
that  it  provides  the  opportunity  to 
view  enhancing  structures  within  a 
high-spatial-resolution  format.  In  this 
format,  morphologic  cues  such  as  stel¬ 
late  patterns  or  irregular  borders  of 
tissues  or  lesions  are  correctly  associ¬ 
ated  with  their  dynamic  enhance- 
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A  Figure  2.  Data  from  a  coronal  image  obtained  through  the  breast  of 
"  56-year-old  patient  with  a  stable  fibroadenoma,  (a)  Early- time-point  (ie, 
precontrast)  keyhole  image,  (b)  Subtraction  image  from  a  late  time  point 
(t  =  355  seconds)  of  dynamic  imaging  shows  contrast  enhancement  in 
the  mass  (KOI),  (c)  Subset  of  time  series  for  this  section  illustrates  a  grad¬ 
ual  contrast  change.  "  =  seconds,  (d)  Time-course  analysis  of  the  mass 
indicates  that  this  fibroadenoma  has  a  moderate  to  long  time  constant 
(t  —  119  seconds). 


merit  properties.  An  alternative  to  the 
keyhole  method  is  rapid  acquisition 
of  low-spatial-resolution  data  sets  that 
are  reconstructed  and  viewed  inde¬ 
pendently  from  high-resolution  data 
sets.  However,  we  believe  the  com¬ 
plexities  of  keyhole  reconstruction  are 
justified  to  ensure  the  temporal  en¬ 
hancement  profiles  and  spatial  prop¬ 
erties  of  a  region  are  correctly  related. 
Spatial  blur  of  enhancement  that  re¬ 
sults  from  acquisition  of  limited  k 
space  will  spread  along  phase-encode 
directions  to  bordering  tissue.  That  is, 
truly  unenhancing  tissues  can  dem¬ 
onstrate  apparent  enhancement  that 
is  actually  caused  by  bleeding  of  sig¬ 
nal  intensity  from  adjacent  enhancing 
tissues.  Vertical  smearing  of  signal 
intensity  from  small  enhancing  ob¬ 
jects  is  illustrated  in  the  subtraction 
images  of  Figures  lb  and  2b. 

In  addition,  the  apparent  amplitude 
of  enhancement  within  truly  en¬ 


hanced  tissue  will  be  diminished.  This 
effect  in  the  phase-encode  dimen¬ 
sions,  y  and  z,  is  demonstrated  with 
the  function  that  characterizes  a  uni¬ 
formly  enhancing  lesion  as 

L(y,z,t)  =  L(y,z,0)  •  [1  +  C(f)] 

=  ^"1{L(ky/kz,0) 

■  [1  +  C(f)]},  (2) 

where  C(t)  defines  the  temporal  prop¬ 
erties  with  C(0)  of  0  and  L(ky,kz,0)  de¬ 
scribes  spatial  properties  of  lesions  in 
k  space.  Changes  in  contrast  in  the 
reconstructed  keyhole  image  of  this 
lesion,  ^(y,z,t),  are  limited  by  the 
spatial  frequencies  that  lie  within  the 
sampled  keyhole  matrix.  Thus  for 
K(kykz)  =  1  within  the  keyhole  range 
and  zero  otherwise,  we  have 

J^(y,z,t)  =  ^9r"1{L(ky,kz,0) 

•  [1  +  K(ky,kz)  •  C(t)]}.  (3) 


The  apparent  enhancement  time 
course  of  the  lesion  is  derived  by  mul¬ 
tiplying  Jz?(y,z,t)  by  an  ROI  mask  and 
then  performing  spatial  integration. 
The  convolution  theorem  allows  us  to 
relate  this  in  k  space  as 

AROI(t)  =  ROI(f)  -  ROI(O) 

=  C(t)  f  R(ky,K) 

■  L(ky,kz) 

•  K^kjdkydk^  (4) 

where  R(ky,kz)  is  the  Fourier  transform 
of  the  ROI  mask.  If  we  assume  the  spa¬ 
tial  pattern  of  the  ROI  closely  matches 
that  of  the  enhancing  lesion,  then 
R(ky/kz)  =  L(ky,kz);  thus,  the  integrand 
becomes  |  L(y,kz)  |2  within  the  keyhole 
area  and  zero  elsewhere.  Finally,  the  full 
matrix  acquisition  is  used  to  normalize 
the  contrast  suppression  effect: 

AROI(f)  s  p-C(t)-  ROI(O); 

where 

P  =  Jxeyholearea  I  k(ky,kz)  | 2  dkydkz 

*  JfuII  k  space  I  L(ky,kz)  | 2  dkydkz .  (5) 

The  factor  (3  represents  the  degree 
of  suppression  in  contrast  change  that 
results  from  sampling  only  a  portion 
of  k  space.  For  the  above  derivation  to 
apply,  the  time  to  complete  each  key¬ 
hole  pass  must  be  short  or  must  be  at 
least  comparable  to  the  period  of 
strong  contrast  change.  The  lesion 
contrast  state  is  considered  nearly 
constant  for  each  keyhole  pass  but  is 
allowed  to  vary  from  pass  to  pass. 
While  coverage  of  full  k  space  may 
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result  in  undersampling  of  temporal 
changes,  such  coverage  is  still  used 
for  normalization  to  provide  a  worst- 
case  bound  for  (3.  Within  the  validity 
of  these  approximations,  sampling  of 
partial  k  space  reduces  the  amplitude 
of  the  enhancement  curve,  although 
the  temporal  shape  of  the  curve  is 
correct.  As  a  result,  time-constant  fit¬ 
ting  of  the  shape  of  the  enhancement 
curve  is  relatively  immune  to  the  ef¬ 
fects  of  partial  sampling  of  k  space  as 
long  as  the  sampling  period  is  ad¬ 
equately  short  and  the  signal-to-noise 
ratio  is  sufficiently  high  to  enable  de¬ 
tection  of  attenuated  changes  in  signal. 

For  illustration,  consider  our  de¬ 
fault  settings,  in  which  one-eighth  of 
k  space  is  acquired  with  the  32  ky  x 
16  kz  keyhole  matrix.  3D  acquisition 
with  a  rectangular  30  x  15-cm  field  of 
view  and  4-mm  sections  gives  phase- 
encode  dimensions  of  1.17  X  4.00-mm 
pixels.  The  change  in  suppression  of 
signal  intensity  for  a  single-pixel  (1.17  x 
4.00  mm)  lesion  would  be  (3  ==  0.13. 
That  is,  if  a  single-pixel  lesion  truly 
changed  contrast  by  100%  (ie,  the  pre¬ 
contrast  level  doubled),  results  with 
this  method  would  indicate  only  a 
13%  change  in  signal  intensity.  For 
larger  lesions,  the  attenuation  effect  is 
lessened  as  the  function  |  L(ky,kz)  |2 
becomes  peaked  within  the  keyhole 
area. 

Examples  of  predictions  with  Equa¬ 
tion  (5)  in  larger  lesions  obtained  with 
the  same  acquisition  parameters  are 
(3  =  0.46  for  a  4.70  x  4.00-mm  lesion, 

P  =  0.81  for  a  8.20  x  8.0-mm  lesion, 
and  p  =  0.92  for  a  11.70  x  12.00-mm 
lesion.  Composite  cinegraphic  display 
of  anatomic  and  subtracted  images 
partially  overcomes  the  impact  of 
contrast  suppression.  Subtraction  re¬ 
moves  time-invariant  image  features 
that  partially  obscure  lesions,  and 
cinegraphic  display  aids  the  observer 
in  distinguishing  true  contrast  changes 
from  noise  and  other  artifact. 

Data  obtained  in  large  lesions  in 
which  p  ~  1  indicate  signal  intensity 
can  be  enhanced  several  times  more 
than  that  depicted  on  nonenhanced 
images,  but  lesions  with  less  change 
in  signal  intensity  are  also  detectable 
and  well  quantified  with  enhance¬ 
ment-time  constants.  Assuming  true 
changes  in  contrast  are  double  the 
precontrast  levels,  we  anticipate  le¬ 
sions  more  than  5  mm  in  diameter 
(where  p  ~  0.5)  would  be  detectable 
with  this  method  since  the  apparent 
enhancement  would  be  -50%  or 
more.  Of  course,  weaker  enhance¬ 
ment  or  increased  background  noise 
or  artifact  would  decrease  conspicuity 
of  these  smaller  lesions. 


The  high  sensitivity  of  contrast- 
enhanced  MR  imaging  for  detection 
of  breast  disease  is  well  known.  Un¬ 
fortunately,  the  modality  has  moder¬ 
ate  specificity  and  is  susceptible  to 
false-positive  findings,  since  some 
benign  tumors  and  tissues  also  exhibit 
enhancement.  Moreover  if  greater 
attention  is  to  be  paid  to  rapidly  en¬ 
hancing  regions,  normal  structures 
need  to  be  appropriately  identified 
that  naturally  exhibit  rapid  signal 
change.  Cinegraphic  images  viewed 
as  individual  sections  (when  acquired 
as  such)  or  drawn  from  a  3D  volume 
set  display  numerous  small  rapidly 
enhancing  regions  representing  blood 
vessels.  To  help  identification  of  ves¬ 
sels  that  can  be  potentially  mistaken 
for  small  lesions,  we  further  exploit 
the  3D  nature  of  the  data  and  pro¬ 
duce  maximum-intensity-proj  ection 
images  of  subtraction  images  obtained 
through  the  full  imaged  volume  along 
orthogonal  directions.  Section-to-sec- 
tion  continuity  of  vessels  is  more  eas¬ 
ily  appreciated  in  this  display.  These 
subtraction  maximum-intensity-pro¬ 
jection  images  are  viewed  cinegraphi- 
cally,  thus  they  also  provide  a  quick 
survey  of  the  full  study.  In  this  format, 
changes  in  spatial  and  temporal  con¬ 
trast  of  a  mass  relative  to  other  breast 
tissues,  vessels,  and  blood  in  the  heart 
(usually  partially  included  in  the  im¬ 
aged  volume)  are  readily  apparent. 
Excerpts  of  a  time  series  of  superior- 
to-inf  erior  maximum-intensity-proj  ec¬ 
tion  images  of  subtraction  images  are 
shown  in  Figure  4a  with  enhance¬ 
ment-time-course  curves  for  regions 
shown  graphically  in  Figure  4b. 

There  are  known  fundamental  tis¬ 
sue  changes  that  provide  the  ratio¬ 
nale  for  the  hypothesis  that  malignant 
breast  cancers  will  preferentially  dis¬ 
play  rapid  changes  in  contrast  over 
most  benign  and  normal  tissues.  Con¬ 
siderable  experimental  data  show  that 
tumors  require  an  enhanced  blood 
supply  to  support  their  rapid  growth 
(26-28).  Thermal  diffusion-based  de¬ 
livery  of  nutrients  sustains  tumor 
nodules  up  to  1-2  mm  in  largest  di¬ 
ameter,  but  their  continued  linear  and 
rapid  exponential  growth  to  larger 
diameters  is  facilitated  by  neovascu¬ 
larity.  Experimental  evidence  that 
tumors  are  angiogenesis  dependent  is 
summarized  in  an  article  by  Folkman 
(29),  in  which  several  studies  are  cited 
that  demonstrate  tumor  size  is  linear, 
slow,  or  stagnant  before  vasculariza¬ 
tion  but  switches  to  rapid  growth  af¬ 
ter  vascularization.  Angiogenesis  does 
not  necessarily  imply  malignancy, 
however,  since  some  benign  tumors 
are  also  angiogenic.  Significant  corre- 
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Figure  3.  Scatter  plot  of  time  constants  de¬ 
rived  from  time-course  data  for  mono-expo¬ 
nential  fits  of  enhancement.  ROIs  were  de¬ 
fined  on  various  tissue  across  12  patients. 
The  N  values  in  the  tissue  type  axis  refer  to 
the  number  of  ROIs  associated  with  each  tis¬ 
sue  type;  only  one  ROI  was  defined  for  each 
discrete  mass.  One  patient  had  two  invasive 
ductal  carcinoma  tumors,  and  another  pa¬ 
tient  had  two  invasive  lobular  carcinoma  tu¬ 
mors.  1  =  invasive  ductal  carcinoma  (n  ~  5), 
2  =  metastatic  neuroendocrine  carcinoma 
(n  =  1),  3  =  invasive  lobular  carcinoma  (n  — 
3),  4  =  enhancing  tissue  adjacent  to  invasive 
ductal  carcinoma  ( n  =  5),  5  =  fibroadenoma 
(n  =  1),  6  =  benign  mass  (n  -  3),  7  =  areolar 
or  nipple  tissue  (n  =  8),  8  =  enhancing  nor¬ 
mal  tissue  (n  =  7).  No  more  than  one  normal 
tissue  ROI  was  derived  from  any  given  patient. 


lations  between  tumor  angiogenesis 
and  invasiveness  and  metastatic  po¬ 
tential  of  breast  cancer  have  been 
demonstrated  with  direct  optical  mi¬ 
crovessel  counts  in  breast  tumor 
specimens  (26-28).  There  is  also  some 
evidence  that  precancerous  breast 
tissues  have  a  greater  angiogenic  pro¬ 
pensity,  suggesting  that  angiogenesis 
has  potential  as  a  precancer  marker 
(30).  Factors  that  trigger  and  moder¬ 
ate  angiogenesis  are  not  fully  known, 
but  their  determination  and  potential 
exploitation  for  diagnosis  and  treat¬ 
ment  are  clearly  critical  areas  for  in¬ 
vestigation. 

Evaluation  of  the  clinical  value  of 
imaging  modalities  sensitive  to  breast 
vascularity,  such  as  MR  imaging  and 
Doppler  US  (31),  is  also  an  active  area 
of  investigation.  The  value  of  tempo¬ 
ral  parameters  derived  from  dynamic 
contrast-enhanced  MR  imaging  in 
helping  distinction  of  benign  and  ma¬ 
lignant  masses  is  an  area  of  debate. 
Currently,  there  is  disagreement 
among  investigators  who  have  per¬ 
formed  clinical  studies  to  evaluate  the 
optimal  method  for  performing  con¬ 
trast-enhanced  MR  imaging  for  the 
characterization  of  breast  masses.  As 
mentioned,  technical  factors  related 
to  volume  of  tissue  coverage,  tempo¬ 
ral  sampling  rate,  and  quantification 
of  temporal  parameters  are  variable 


140  •  Radiology 


July  1995 


Sample  Enhancement  Curves 


Time  (sec) 


a.  b- 

Figure  4.  (a)  Maximum  intensity  projection  of  subtraction  images  in  the  superior-inferior  projection.  Six  noncontiguous  time  frames  ”e 

shown  from  images  normally  viewed  cinegraphically.  ROI#3  delineates  an  invasive  ductal  carcinoma.  Enhancement  in  linear  structures  delin- 
ea°«ves”  = anterior?  L  =  left,  P  =  posterior,  R  =  right,  f  =  time  (seconds).  "  =  seconds,  (b)  Corresponding  enhancement-hme-course 
curves  for  ROIs  illustrate  changes  in  relative  rate  and  intensity  of  enhancement.  Curves  were  derived  from  ROIs  drawn  on  related  single- 
section  cinegraphic  displays.  ROI#l  =  left  ventricular  blood  (broken  line),  ROI#2  =  liver  (O),  ROJ#3  -  invasive  ductal  carcinoma  (solid  lin  ), 
ROI#4  =  area  adjacent  to  mass  (dashed  line). 


across  these  studies.  The  immediate 
technical  objectives  of  this  study  were 
to  overcome  these  critical  limitations. 
These  objectives  are  largely  met  with 
the  3D  keyhole  technique,  which  of¬ 
fers  full  tissue  coverage  at  high  tem¬ 
poral  sampling  for  subsequent  quantita¬ 
tive  enhancement  time  characterization. 
As  we  have  discovered,  the  ability  to 
dynamically  study  entire  breasts  at  MR 
imaging,  thus  obviating  prospective 
identification  of  one  or  all  lesions  amidst 
complex  contrasts,  is  a  great  practical 
advantage.  A  drawback  of  this  approach 
is  the  large  data  volume  generated  per 
case.  Currently,  we  resolve  the  data 
management  issue  by  performing  recon¬ 
struction,  display,  and  analysis  on  an 
off-line  workstation  with  customized 
software  tools.  Our  software  has  been 
specifically  designed  to  present  data  in  a 
variety  of  formats  for  specific  purposes, 
such  as  the  survey  of  all  tissues  in  cine- 
graphic  maximum-intensity-projection 
images  of  subtraction  images  or  the  cine- 
graphic  viewing  of  a  single  section  for 
quantitative  analysis  of  individual  tis¬ 
sues.  This  method  also  benefits  from  the 
other  well-known  features  of  breast  MR 
imaging,  such  as  the  ability  to  image 
breasts  with  dense  parenchyma  or  im¬ 
plants  and  to  provide  high-quality  refor¬ 
matted  tomographic  views  that  can  po¬ 
tentially  aid  planning  for  surgery.  ■ 
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